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We study the impact of neutrino oscillations on the interpretation of the supernova (SN) 1987A 
neutrino signal by means of a maximum-likelihood analysis. We focus on oscillations between F e 
with 77 M or 77 T with those mixing parameters that would solve the solar neutrino problem. For the 
small-angle MSW solution (Am 2 « 1CT 5 eV 2 , sin 2 29o ~ 0.007), there are no significant oscillation 
effects on the Kelvin-Helmholtz cooling signal; we confirm previous best-fit values for the neutron- 
star binding energy and average spectral V e temperature. There is only marginal overlap between 
the upper end of the 95.4% CL inferred range of (Ep e ) and the lower end of the range of theoretical 
predictions. Any admixture of the stiffer F M spectrum by oscillations aggravates the conflict between 
experimentally inferred and theoretically predicted spectral properties. For mixing parameters in the 



neighborhood of the large-angle MSW solution (Ar 



10 _5 eV 2 , sin 2 



0.7) the oscillations 



in the SN are adiabatic, but one needs to include the regeneration effect in the Earth which causes 
the Kamiokande and 1MB detectors to observe different V e spectra. For the solar vacuum solution 
(Am 2 « 10" 10 eV 2 , sin 2 29 « 1) the oscillations in the SN are nonadiabatic; vacuum oscillations 
take place between the SN and the detector. If either of the large-angle solutions were borne out by 
the upcoming round of solar neutrino experiments, one would have to conclude that the SN 1987A 77 M 
and/or V e spectra had been much softer than predicted by current treatments of neutrino transport. 
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I. INTRODUCTION 

Neutrino oscillations can modify the characteristics of 
the neutrino signal from a supernova (SN), in particular 
if matter effects are included |5J]. After the observation 
of the SN 1987A neutrinos by the Kamiokande || and 
1MB H detectors many authors Q] discussed the impact 
of matter-induced oscillations on the prompt v e burst be- 
cause the first event at Kamiokande had been observed 
in the forward direction, allowing for an interpretation in 
terms of v e -e scattering. If this interpretation were cor- 
rect one could exclude a large area of mixing parameters 
where the MSW effect in the SN envelope would have 
rendered the prompt v e burst unobservable. 

Because a single event does not carry much statistically 
significant information (the first Kamiokande event may 
have coincidcntally pointed in the forward direction), a 
more interesting question for the interpretation of the 
SN 1987A neutrino signal is the impact of oscillations 
on the main v e pulse which is detected by the reaction 
V e p — » ne + . The SN emits roughly equal amounts of 
energy in (anti)neutrinos of all flavors, but with different 
spectral characteristics. Current treatments of neutrino 
transport yield 
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i.e. (E Vs ) w l(Ev e ) and {E v ) « |(£^ e ) for the other fla- 
vors. A partial conversion between, say, F^'s and V e 's due 
to oscillations would "stiffen" the V e spectrum observable 
at Earth J8|,|9| . (We will always take V e -V li oscillations to 
represent either Ve-v^ or v e -v r oscillations.) Within a 
plausible range of progenitor star masses and depending 
on the equation of state, numerical computations yield 



E h = 1.5 -4.5 x 10 53 erg 



(2) 



for the total amount of binding energy . It is almost 
entirely released in the form of neutrinos. 

The expected average SN 1987A V e energy implied by 
the detected signal is about 9— lOMeV, with a 95.4% 
confidence interval reaching up to 14 MeV in some anal- 
yses [ ^0|JTl| |, i.e. barely up to the lower end of the the- 
oretical predictions quoted in Eq. (Q). If a partial swap 
v e «-> Vy, had occurred, the expected F e energies should 
have been lower, causing an even larger strain between 
measured and predicted j/ e energies. For an "inverted" 
mass matrix with m Ue > m„ the V e -V^ oscillations would 
have been resonant and thus nearly complete for a large 
range of mixing parameters. Therefore, such inverted- 
mass schemes are likely excluded on the basis of the 
SN 1987A data [|jl|. 

If the mass hierarchy is "normal" with m„ c < m v , 
oscillations in the antineutrino sector are significant only 
for large mixing angles which are often thought to be un- 
likely. Therefore, in the original analyses of the SN 1987A 
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neutrinos, little attention has been paid to antineutrino 
oscillations. 

Since then much progress has been made with the ob- 
servation of solar neutrinos in four experiments which all 
report a deficit and thus point to oscillations. While it 
remains uncertain if the solar neutrino deficits are indeed 
caused by oscillations, it has become clear that there is 
no simple "astrophysical solution." If the oscillation in- 
terpretation is adopted there remain three islands in the 
sin 2 20o-Am 2 -plane (vacuum mixing angle Oo) where 
the results from all experimental measurements of the so- 
lar neutrino flux are consistently explained, namely the 
"vacuum solution" with Am 2 near 10 _10 eV 2 and nearly 
maximum mixing the "small-angle MSW solution" 
with Am 2 around lCT 5 eV 2 and sin 2 26 « 0.007, and 
the "large-angle MSW solution" with about the same 
Am 2 and sin 2 20o in the neighborhood of 0.7 S. It will 
turn out that if one of the large-angle solutions would 
be borne out by one of the forthcoming experiments Su- 
perkamiokande, SNO, or BOREXINO, then a significant 
impact on the interpretation of the SN 1987A signal could 
not be avoided. 

In a recent study, Smirnov, Spergel, and Bahcall ||] 
found that the large-angle solutions were essentially ex- 
cluded by the SN 1987A data because of the "stiffened" 
spectra they would have caused at the detectors. How- 
ever, this conclusion relies heavily on theoretical predic- 
tions for the spectral properties of a SN neutrino signal. 
Kernan and Krauss ]14|] , on the other hand, arrive at 
the opposite conclusion, namely that a significant oscil- 
lation effect was actually favored by the data. Of course, 
they discard certain theoretical predictions for the sig- 
nal characteristics. Smirnov, Spergel, and Bahcall have 
performed a joint analysis for the Kamiokande and 1MB 
detectors. However, in the neighborhood of the large- 
angle MSW solution, matter-induced oscillations in the 
Earth are important. They cause a different amount of 
"regeneration" of the oscillations on the neutrino path 
through the Earth which was 3900 and 8400 km for the 
Kamiokande and 1MB detectors, respectively, which thus 
would have observed different v e spectra [^5[. Kernan 
and Krauss, on the other hand, have only considered 
nonadiabatic oscillations which restrict the validity of 
their analysis to Am 2 < 10~ 10 eV 2 , thus ignoring the 
important case of the large-angle MSW solution. 

Therefore, we presently reexamine the impact of large- 
angle neutrino oscillations on the SN 1987A signal in- 
terpretation. If neutrino oscillations between v e and an- 
other flavor occur at all with a large mixing angle, the 
mixing parameters probably correspond to those solving 
the solar neutrino problem. Therefore, we focus on mix- 
ing parameters in the neighborhood of the large-angle 
MSW solution and of the vacuum solution of the solar 
neutrino problem. We will assume thermal neutrino spec- 
tra with different temperatures for the v^s and F e 's. We 
will then perform a maximum-likelihood analysis for the 
neutrino temperature and total emitted energy. 

In Sect. II we discuss the assumed primary neutrino 



spectra and their modification by oscillations. Sect. Ill 
is devoted to our statistical methodology and Sect. IV to 
detailed numerical results. In Sect. V we summarize our 
findings. 



II. NEUTRINO SPECTRA 
A. Primary Spectra 

The most detailed statistical analysis of the SN 1987A 
neutrino signal has been performed in the papers by 
Loredo and Lamb where one of the main goals 

was to estimate the Kelvin-Hclmholtz cooling time scale 
of the newly formed neutron star, and to derive limits on 
the V e mass from the absence of pulse dispersion effects. 
Therefore, the time structure of the neutrino signal was 
crucial; it had to be parametrized in terms of a variety 
of cooling models. In our study, on the other hand, we 
will focus on the spectral characteristics of the neutrino 
fluence (time-integrated flux) and their modification by 
oscillations. Because we will need to vary neutrino mass 
differences and mixing angles, the overall number of pa- 
rameters would get out of hand if we were to analyse 
the time structure of the burst together with neutrino 
oscillation effects. 

Numerical simulations ]rj| and an analytic argument 
|]l9f indicate an approximate equipartition of the energy 
emitted in different (anti)neutrino species with different 
time-averaged energies as quoted in Eq. (Q) . The detailed 
spectral shape, however, is not well known. Monte-Carlo 
studies of neutrino transport indicate that the instan- 
taneous neutrino spectra are "pinched," meaning that 
their low- and high-energy parts are suppressed relative 
to a Maxwcll-Boltzmann spectrum of the same average 
energy. Usually the instantaneous spectra are expressed 
in the form |16|] 

where 77 is an effective degeneracy parameter. Both T 
and r\ are functions of time. It must be stressed that 
the "pseudo degeneracy parameter" 77 for and v T is 
the same as that for and 77 r , in contrast with the 
degeneracy parameter of a real Fermi-Dirac distribution 
which has the opposite sign for antineutrinos relative 
to neutrinos. Therefore, Eq. (^|) is a somewhat arbi- 
trary two-parameter representation of the neutrino spec- 
tra which allows one to fit two of their moments, for 
example (E) and (E 2 ). Janka and Hillebrandt |l6| found 
that throughout the emission process 77 decreases from 
about 5 to 3 for v e , from about 2.5 to 2 for F e , and from 
about 2 to for u^, T and v^^- 

The time-integrated spectrum, however, need not be 
pinched. We characterize it by the moments (E) and 
(E 2 ), and call it "pinched" if the ratio (E 2 )/(E) 2 is 
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smaller than for the Maxwell-Boltzmann case, "anti- 
pinched" otherwise. As a simple example we consider 
a cooling model where neutrinos are emitted from a neu- 
trino sphere with a fixed radius and an exponentially 
decreasing effective temperature. If the instantaneous 
spectra are of the form Eq. (||) with a fixed rj, then the 
time-integrated spectrum is pinched for 77 ^ 1.7 and an- 
tipinched for rj < 1.7. For rj w 1.7 it is approximately of 
the Maxwell-Boltzmann form. 

An exponential cooling model is, of course, very sim- 
plistic. In a real SN the v e temperature will initially 
rise, and may stay approximately constant for some time, 
while the effectively radiating surface shrinks quickly 
within the first second. Still, the exponential cooling 
example illustrates that a thermal Maxwell-Boltzmann 
spectrum may be a relatively good approximation for the 
time-integrated spectrum because of the compensating 
effects between instantaneous pinching and the superpo- 
sition of different spectra in the course of the protoneu- 
tron star's cooling history. Certainly, there is no reason 
to expect the time-integrated spectrum to be of the form 
Eq. (03). This parametrization does not allow one to de- 
scribe antipinched spectra, only pinched ones. 

For the rest of this study we will make the simplify- 
ing assumption that the time-integrated spectra are de- 
scribed by the Maxwell-Boltzmann form 



F{E) 



dtf{E,t) oc E 



2-E/T 



(4) 



with a different effective temperature for T> e and v^. 
These "temperatures" are parameters which characterize 
the time- integrated spectra by virtue of T = \{E) and 
thus do not exactly correspond to a physical temperature 
at the neutron star. 



B. Modification by Oscillations 

In the Kamiokande and 1MB detectors, SN neutrinos 
are almost exclusively detected by the reaction v e p — > 
ne + where the final-state positron is measured by its 
Cherenkov emission of photons. If neutrinos do not mix, 
their fluence F„ c (E) relevant for the detection process is 
identical with the pimary v e spectrum F^ (E) emitted 
from the SN. In the presence of v e <-» oscillations, on 
the other hand, each primary 77^ arrives with a proba- 
bility p in the V e flavor state at the detector, while each 
primary v e arrives as V e with the "survival probability" 
1 — p so that 



F Vc = (l-p)Fl+pF» 



(5) 



This incoherent superposition of the individual flavor 
fluxes is justified by the incoherent neutrino emission 
from different regions in the star and by different pro- 
cesses |J. 

The "permutation factor" p is in general a function of 
the neutrino energy E, the mass difference Am 2 , and the 



vacuum mixing angle Qq. In addition, it is important to 
note that the neutrinos are produced in a region of high 
matter density. The effective mixing angle in a medium 
is given by the well-known formula 



tan 29 



sin 26n 



cos29o T p/p 



(6) 



where 0o is the vacuum mixing angle, p the matter den- 
sity, and the upper sign refers to v, the lower to V. The 
"resonance density" is defined by 



mjvAm 2 
2^2 G F Y e E' 



(7) 



where Am 2 = 

mixture of and m\ that of v, 



2 — m\ with mi the dominant mass ad- 
For neutrinos with 



a normal mass hierarchy (mi > mi) the denominator 
in Eq. (Q) vanishes for p = p rcs cos29 , causing max- 
imum mixing with O = 7r/4 and thus a "resonance." 
For antineutrinos, and because we always assume a nor- 
mal mass hierarchy, the denominator of Eq. (|^) is always 
larger than cos 20o so that the medium mixing angle is 
always smaller than the vacuum one. 

For our purposes with neutrino energies E > lOMeV 
and mass differences Am 2 < 1CP 3 eV 2 the resonance den- 
sity is of order 10 3 gem -3 or less. With p w 10 12 gem -3 
at the neutrino sphere, the effective antineutrino mixing 
angle at the production site is < 10~ 9 , even if the vac- 
uum mixing angle is maximal. Therefore, the medium ef- 
fects "demix" the antineutrinos, causing the flavor eigen- 
states at the production site to coincide essentially with 
the propagation eigenstates. 

As the neutrinos leave the SN they propagate through 
a certain density profile and ultimately reach the sur- 
rounding vacuum. The Am 2 values corresponding to the 
large-angle solutions of the solar neutrino problem are 
representative of two cases that need to be distinguished 
for the further flavor evolution of the neutrino burst. 

The simpler case is the vacuum solution for Am 2 < 
10~ 10 eV 2 . The propagation out of the SN is not adia- 
batic so that the neutrinos emerge essentially as flavor 
eigenstates which then oscillate on their way to Earth. 
Therefore, the permutation factor has the form 



P 



\ sin 2 26 . 



(8) 



We note that Am 2 w 10~ 10 eV 2 is at the borderline for 
this statement to apply; for slightly larger mass differ- 
ences the detailed propagation through the SN envelope 
must be taken into account || . 

For the large-angle solar MSW solution with Am 2 w 
10 -5 eV 2 we are in the adiabatic regime where the neu- 
trinos stay in a propagation cigenstate throughout their 
journey out of the SN ||. What emerges is a flux of rti\ 
eigenstate neutrinos with the v e spectrum, and one of mi 
eigenstates with the spectrum. 

We stress that this statement applies even though the 
neutrinos encounter a density discontinuity correspond- 
ing to the outward moving shock wave which ultimately 
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ejects the SN mantle and envelope. At the neutrino 
sphere, the propagation and flavor eigenstates coincide 
because of the medium-induced demixing effect described 
above. When the neutrinos encounter a density discon- 
tinuity in a medium so dense that they are sufficiently 
demixed, then no significant flavor transitions will occur 
even though this discontinuity violates the adiabaticity 
condition. Within the first few seconds after collapse the 
shock wave may reach a radius of at most a few 10 5 km. 
In typical progenitor star models the density varies ap- 
proximately as r -3 . Initially, the neutrino sphere with 
a density of about 10 12 gcm -3 is at a radius of about 
100 km. Therefore, within the Kelvin-Hclmholtz cooling 
phase the shock wave may reach a density about 9 or- 
ders of magnitude smaller than the neutrino sphere, i.e. 
a density as low as 10 3 gcm~ 3 . For Am 2 w 10 -5 eV 2 the 
resonance density is about 10 gem -3 . Hence, during the 
entire Kelvin- Helmholtz cooling phase the medium mix- 
ing angle is small when the neutrinos encounter the shock 
wave. Therefore, the impact of level crossing between the 
propagation eigenstates on the neutrino spectra arriving 
at the detector can be neglected. 

Because neutrinos with Am 2 w 10 -5 eV 2 emerge from 
the SN as propagation eigenstates, no oscillations occur 
on the way from the SN to Earth. Thus, we would have 
p — sin 2 0o if there were no further intervening matter. 

However, in order to reach the Kamiokande and 1MB 
detectors, the neutrinos had to traverse c?kam = 3900 km 
and diMB = 8400 km of matter in the Earth, with an 
average density of about £kam = 3.4 gem -3 and pimb = 
4.6 gem -3 , respectively Therefore, the permutation 
factor relevant for each detector is || 

p = sin 2 9 - sin 29 sin(26 - 29) sin 2 (nd/£) . (9) 

The medium mixing angle relevant for each detector is 
given by Eq. (||) with p = p K AM or pi M B, respectively, 
the distance in Earth is d = g?kam or c?imb, and the 
oscillation length is 



£ = 



4nE sin 29 



Am 2 sin29 



(10) 



with the relevant medium mixing angle. For the solar 



vacuum solution with Ar 
is unimportant. 
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III. STATISTICAL METHODOLOGY 



A. Parameter Estimation and Confidence Regions 



which exploits the sparse data efficiently. The maximum- 
likelihood method [2^,||l[ is particularly well suited for 
such problems, i.e. problems where it is essential to ex- 
tract the maximum possible information from a small 
number of events. This method has been used by several 
authors to analyse the SN 1987A neutrino signal, e.g. 
Refs. |0|Jll]Jl|Jl|]. 

The method consists of deriving the set of parameters, 
collectively denoted by a, for which the probability of 
producing the observed data set, collectively denoted by 
x, becomes maximal. The probability density as a func- 
tion of a for producing the observed data is called the 
likelihood function C(x, a). The maximum-likelihood es- 
timation a* for the true but unknown parameter set «o 
is implicitly defined by 



C(x, a*) = max C(x, a), 



(11) 



where D is the parameter domain. 

An estimation a* of the true parameters ao is useful 
only if one also determines a confidence region around 
a* which contains the true parameters with a specified 
probability (3. To construct this region assume that the 
true parameters «o are given. We can then determine 
the probability distribution P Qo (a*) of the likelihood es- 
timator and define a region Dp_ ao from the condition 
P<x ( a *) > P f° r a * G Df3 Ma . To make it unique we ad- 
ditionally require that P ao (a*) is larger for all a* within 
Dp,a than for those outside. Put another way, we re- 
quire -D/3, Qo to be bounded by a contour of constant 
Pa ( a *)- The confidence region can now be defined 
as the region of parameters a for which a* € Dp a . Note 
that this set is in general not equal to Dp^ air . 

In practice, this region is difficult to calculate because 
finding Dp a alone requires integrating over the space of 
possible observations, a task usually achieved by Monte- 
Carlo sampling. However, if C is Gaussian the confidence 
region is given by the condition 



ln£(a;, a*) — \nC{x, a) < \ Xp(k), 



(12) 



again with the additional requirement that it should be 
bounded by a contour of constant C in parameter space 
pof . Further, k is the number of parameters which for our 
study will usually be k = 2. Note that x (2) = 2.3, 4.61, 
and 6.17 for /3 = 68.3%, 90%, and 95.4%, respectively. 
We stress that the confidence regions thus determined 
are not exact, especially when they are very distorted so 
that the parameters are strongly correlated. 



The purpose of the present study is to estimate the 
parameters E\> and Tp e which characterize the neutrino 
fluence from SN 1987A and to study the impact of neu- 
trino mixing on this estimate. Because of the small num- 
ber of SN 1987A events in the Kamiokande and 1MB 
detectors this task is rather delicate. One needs a sta- 
tistical estimator which is consistent and unbiased, and 



B. Likelihood Function 

It is not trivial to determine the likelihood function ap- 
propriate for our problem. The primary observations of 
the water Cherenkov detectors consist of the information 
when a given photomultiplicr has fired. This information 
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can be used to reconstruct the event location in the de- 
tector and the energy of the detected charged particle. 
For our purposes it is probably sufficient to use the re- 
ported event energies as the primary data set and assume 
that they are related to the true positron energies by a 
Gaussian distribution. 

In order to model the likelihood function we consider 
detection energy bins [Ei,Ei + SE] with i = 1,, .., Nun- 
The spectrum of detected energies is n(E) so that the 
number of expected counts in bin i is to lowest order 
n(Ei) 5E. However, in a real experiment one obtains 
an integer number 2Vj of counts in a given bin i. The 
probability for such an outcome is 



_ [njE^SE]** (Ei)SE 



NA 



(13) 



where the Nt are the actual observations and thus repre- 
sent the data. The likelihood function is 



C 



N bin 
i=l 



This expression can be transformed to 



N oba 



C = d 



(14) 



(15) 



where N \^ a is the total number of experimentally ob- 
served events. The constant C is irrelevant for the 
purpose of parameter estimation and the determina- 
tion of confidence regions. For a joint analysis of the 
Kamiokande and 1MB detectors, the likelihood function 
is the product of the likelihood functions for each detec- 
tor. 



C. Expected Energy Spectrum 

In order to translate the V e fluence F^ e (E) at Earth 
to an expected spectrum n(E) of counts we must first 
determine the energy spectrum of secondary positrons 
in the V e p — > ne + reaction. Its cross section Ov fV as a 
function of neutrino energy E is 



0"O 



E 

in,. 



'"I 



2Q , Q 2 -m 2 



E 



E 2 



1/2 



(16) 



where Q = 1.29 MeV is the neutron-proton mass differ- 
ence, m e the electron mass, and cto = 2.295 x 10 _44 cm 2 . 
We ignore Coulomb and radiative corrections as well as 
neutron recoils. Therefore, the positron spectrum in the 
detector is 



n+(E) 



4ttD 2 



(E + Q)F Ve {E + Q), 



(17) 



where D — 50 kpc is the distance to the SN and N p the 
number of target protons in a given detector, namely 
1.43 x 10 32 for Kamiokande and 4.55 x 10 32 for 1MB. 



The positron spectrum n + (E) produced in the detec- 
tor is not identical with the spectrum n(E) of events that 
one expects to detect. The reported energy -Edet for an 
event is reconstructed from the number of photomultipli- 
ers that have been triggered by the Cherenkov light of the 
positrons produced in the detector. Because this involves 
a Poissonian process, a certain number of active photo- 
multipliers corresponds to a range of possible positron 
energies E + that may have caused this event. Moreover, 
there is an E + dependent efficiency curve r] (E + ) that a 
given positron will trigger the detector at all. While this 
function is essentially a step function for the Kamiokande 
detector, it is fairly nontrivial for 1MB where about a 
quarter of the photomultipliers were not operational at 
the time of SN 1987A due to a failed power supply. 

The spectrum of possible reconstructed event energies 
Edct that may be attributed to a true positron energy 
E + is not universal throughout the detector; there are 
nontrivial geometry effects. Still, we use a universal dis- 
tribution for the probability of finding E dct if the true 
energy was E + , 



P(E dct ,E+) 



2-Ka{E^ 



exp 



Edct) 2 



2a 2 (£ + ) 



(18) 



Motivated by the Poissonian nature of the detection pro- 
cess we approximate the energy- dependent dispersion by 



a(E+) = \J E a E^ 



(19) 



For each detector we fit E a from the uncertainties of the 
reported experimental event energies We find that 

a good approximation is E„ = 0.75 MeV for Kamiokande 
and 1.35 MeV for 1MB. 

Instead of using a universal function for P(E dct , E + ) 
we could have used the reported experimental errors Oi 
for each event. This procedure would leave our results al- 
most unchanged while causing complications for the def- 
inition of an overall detector efficiency curve below. 

In both detectors a trigger threshold for the minimum 
number of photomultipliers was used in order to attribute 
a given event to an external signal rather than to back- 
ground. This corresponds to a lower E dct threshold of 
E cut = 7.5 MeV for Kamiokande and 19 MeV for 1MB. 
The published trigger efficiency curves rj(E + ) are thus to 
be interpreted as 



t 1 (E+) = m {E+) / dE dct P(E det ,E + ), 



(20) 



where i]q(E + ) represents efficiency reductions from other 
causes such as geometry and dead-time effects. 

In Fig. [I] we show rj(E + ) and rj (E + ) for both 
Kamiokande and 1MB where for the latter detector a 
13% dead-time effect is not taken into account in the 
efficiency curve. For Kamiokande, rjo(E + ) is essentially 
constant down to the threshold, revealing that the effi- 
ciency curve r](E + ) is dominated by the trigger threshold 
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FIG. 1. Efficiency curves for Kamiokande and 1MB. A 
13% dead-time effect for 1MB is not included. The rj 
curves (dashed) represent the overall efficiencies published in 
Refs. while the 770 curves (solid) are corrected accord- 

ing to Eq. ( pp| ) for the "smearing-out" of -Bdet relative to the 
positron energy E+. 



IV. NUMERICAL RESULTS 



A. No Mixing 

For comparison with previous work we begin our max- 
imum-likelihood analysis with the case of no neutrino 
mixing. We search for the best-fit SN binding energy 
and the effective V e temperature which charac- 
terizes the assumed Maxwell-Boltzmann v e spectrum of 
the time- integrated flux by virtue of (Ez; c ) = 3T,7 e . We 
assume equipartition of the released SN energy between 
all (anti)neutrino species so that E\, is given by six times 
the inferred total energy emitted in F e 's. 

In Fig. H we show the contours of constant likelihood in 
the Tj7 e -i5 b -plane which correspond to 68.3%, 90%, and 
95.4% confidence regions, respectively, and the best-fit 
values for Tp e and E^. In the upper panel we show the 
results from separate analyses for the Kamiokande and 
1MB detectors, in the lower panel from a joint analy- 
sis. Our best-fit values for the Kamiokande detector are 



and by the Poissonian nature of the detection process. 
For 1MB, on the other hand, there is a significant geo- 
metrical efficiency modification. 

The expected spectrum of detected energies is thus re- 
lated to the actual positron spectrum by 



n(-Edet) 



dE + P(E det ,E+) m (E + )n + (E + ) (21) 



for Ed c t > E cut , and n(Edet) = otherwise. With this 
result we are armed to perform the maximum likelihood 
analysis. 



D. Detector Background 

The statistical analysis described above ignores the de- 
tector background, i.e. the fact that any event ascribed 
to the SN burst can also be due to background, and con- 
versely, any event attributed to background can have 
been caused by the SN burst. In Loredo and Lamb's 
analyses [[lO 11 the background spectrum was included 
in the expected event rate. Events much earlier or much 
later than the main burst are automatically discrimi- 
nated against and thus do not overdominate the low- 
energy part of the expected event distribution. Without 
the possibility to discriminate against background events 
by the temporal relationship to the main burst we must 
use the cut represented by the energy threshold E cut . 
We stress that including the background as in Loredo 
and Lamb's analyses does not cause a large modification 
of the implied SN binding energy and neutrino tempera- 
ture. 
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FIG. 2. Contours of constant likelihood which correspond 
to 68.3%, 90%, and 95.4% confidence regions, and best-fit 
values for Tp- e and E\>. Upper panel: Kamiokande and 1MB 
separately. Lower panel: Joint analysis. Dashed lines mark 
the 68.3% confidence regions of the separate fit. 
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TABLE I. Best-fit values for the SN 1987A parameters for 
three neutrino mixing scenarios with a relative 77 M tempera- 
ture r = T— /T- e — 2.0 each. The expected event numbers 
and energies result from the joint analysis for the Kamiokande 
and 1MB detector. The maximum likelihood A ln(£ max ) is 
relative to the no-mixing case. The case of vacuum oscilla- 
tions corresponds to Am 2 < 10~ 10 eV 2 but is otherwise inde- 
pendent of the mass difference. 



Tjr e = 2.5 MeV and E h = 4.9 x 10 53 erg while for 1MB 
they are 3.7 MeV and 5.4 x 10 53 erg, respectively. With 
the Kamiokande best-fit spectrum we find 11 neutrino 
events for Kamiokande and about 1 for 1MB. Conversely, 
the 1MB best-fit spectrum yields about 24 Kamiokande 
and 8 1MB events. 

While the overlap between the separate confidence 
contours is somewhat marginal, it is sufficient to al- 
low for a joint analysis. The joint best-fit values are 
T Fe = 3.6 MeV and E h = 3.4 x 10 53 erg. These best- 
fit parameters as well as the event numbers and average 
event energies corresponding to them are summarized in 
Tab. | 

Our results differ somewhat from those of Janka and 
Hillebrandt |l6| in that these authors find more restric- 
tive confidence contours. We believe that the difference 
is caused by their use of a simplified likelihood function 
where E^ot is identified with E + without allowing for a 
smearing-out effect, and by their use of a Gaussian rather 
than a Poissonian modulation of the detection process. 

The inferred neutron-star binding energy agrees well 
with theoretical expectations of E\> — 1.5 —4.5 x 10 53 erg. 
The best-fit (Ep s ), however, is rather low compared with 
the range of theoretical predictions quoted in Eq. (Q); 
only the 95.4% confidence region slightly touches the pre- 
dicted range. 



B. Vacuum Oscillations 

Next, we study the case of vacuum oscillations which 
is relevant for small neutrino mass differences (Am 2 < 
1CP 10 eV 2 ). The swap probability p is given by the sim- 
ple formula Eq. (^) which depends only on the vacuum 
mixing angle so that no explicit dependence on Am 2 



obtains. In analogy to the F e 's we describe the time- 
integrated V y, flux by a Maxwell-Boltzmann spectrum 
with the same total energy, but with a higher effective 
temperature Tp — r Tp e , where the factor r is predicted 
to lie in the range 1.4-2.0. 
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FIG. 3. Maximum likelihood, binding energy, and 77 e tem- 
perature as functions of the vacuum mixing angle. The 77 M 
temperature is given by Tp = rT- a with the indicated r 
values. 
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We begin by performing the maximum-likelihood anal- 
ysis for a fixed vacuum mixing angle and a fixed r-factor 
while allowing E^ and T„ B to float. In Fig. || we show as 
a function of sin 2 20o the maximum likelihood and the 
best- fit Eb and T& e . We show these curves for r = 1.4, 
1.7, and 2.0. 

For r = 2.0 our results agree well with those of Ker- 
nan and Krauss 14 . The maximum-likelihood curve has 
a maximum for sin 20o sa 0.5 so that a relatively large 
mixing angle appears to be favored by the data. The in- 
ferred SN parameters and expected detector signals for 
this case are summarized in Tab. |. In Fig. || the inferred 
best-fit binding energy is greater for large mixing angles 
compared to the no-mixing case, while the best-fit spec- 
tral temperature is a monotonically decreasing function 
of sin 2 26 . For sin 2 26 > 0.5 and r = 2.0 the best-fit 
(EjjJ is below 6MeV. Such a value is far below what is 
predicted theoretically so that it looks like large mixing 
angles are difficult to reconcile with the SN 1987A data. 

We can also fix the binding energy and neutrino tem- 
perature according to theoretical predictions. Figure |] 
shows ln(£) for E h = 3 x 10 53 erg and {Ej, e ) = 14MeV 
as a function of the mixing angle for several values of 
the relative temperature. The likelihood is a mono- 
tonically decreasing function of sin 2 29o so that, tak- 
ing the predicted SN parameters seriously, the best-fit 
mixing angle is zero, and large mixing angles are dis- 
favored. For r = 1.4 the 95.4% confidence interval is 
< sin 2 26 < 0.17. 




0.2 0.4 0.6 0.8 
sin 2 20 o 

FIG. 4. Likelihood for a fixed E h = 3 x 10 53 erg and 
(E- c ) = 14MeV as a function of the vacuum mixing angle. 
The temperature is given by T-p = rT- c with the indi- 
cated r values. 

Suppose that future experiments will establish vacuum 
oscillations as a solution of the solar neutrino problem. 
What would this imply for the SN 1987A parameters? 
To study this question we show in Fig. || the 95.4% con- 
fidence contours in the T^-Sb-plane for a joint analy- 
sis between the detectors with sin 2 20o = 1 and with 
r = 1.0, 1.4, 1.7, and 2.0. 
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FIG. 5. Best-fit values for T^ e and Eh, and contours of 
constant likelihood which correspond to 95.4% confidence re- 
gions. In each case a joint analysis between both detectors 
was performed with sin 2 2Qq = 1 and the indicated relative 
77 M temperature r. The hatched region corresponds to the 
theoretical predictions of Eqs. (Q) and (^|). 

The 1.0 case corresponds to no mixing; the contour 
is identical with that of the lower panel of Fig. ||. The 
maximum v e temperature within the 95.4% confidence 
region is about 4.6 MeV, corresponding roughly to the 
lower limit for the range of predicted (E„ e ) values as 
given in Eq. (|j). 

For t = 1.4 the 95.4% CL region for the v e energies 
does not overlap with theoretical predictions. Therefore, 
if the vacuum solution would be borne out by future solar 
neutrino experiments, one would be forced to conclude 
that there is a significant problem with the predicted SN 
neutrino spectra and energies. 



C. Adiabatic Oscillations and Earth Effect 



The most complicated case obtains if the solar neutrino 
problem is solved by large-angle MSW oscillations where 
Am 2 k, 10 -5 eV 2 . The propagation out of the SN is adia- 
batic so that no oscillations occur between there and the 
Earth, but we need to include regeneration effects caused 
by the matter-induced oscillations in the Earth. The per- 
mutation factor Eq. (||) is different for the two detectors; 
it is a function of the mass difference, the vacuum mixing 
angle and the neutrino energy. 

As in Sect. IV. B we begin by performing the maxi- 
mum-likelihood analysis for a fixed Am 2 and sin 2 20o 
while allowing E\> and T-y e to float. In Fig. ^ we show 
contours of ln(£ max ) relative to the no-mixing value 
ln(£ max ) = —41.0 in steps of 1. We have used flu- 
ences with the same total energy as for z7 e and a rel- 
ative temperature r = 2.0 (upper panel) and r = 1.4 
(lower panel). The shaded areas correspond to a negative 
Aln(£ max ) and thus to a reduced likelihood relative to 
the no-mixing case. We emphasize that these areas can- 
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FIG. 6. Contours of A]n(£ m ax), which is the maximum 
likelihood relative to the no-mixing value ln(£ max ) = —41.0. 
The contour lines are in steps of 1. Shaded areas correspond 
to Aln(£ max ) < 0, i.e. regions which are disfavored relative 
to the no-mixing case. The relative 77 M temperature r was 2.0 
(upper panel) and 1.4 (lower panel). 
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FIG. 7. Contours of Aln(£ m ax) relative to the no-mixing 
case for a fixed SN binding energy Eh — 3 x 10 53 erg. The 
contour lines are in steps of 1. Shaded areas correspond to 
A ln(£ max ) < 0, i.e. regions which are disfavored relative to 
the no-mixing case. The relative z7 M temperature r was 2.0 
(upper panel) and 1.4 (lower panel). 



not be interpreted as being excluded even though they 
are disfavored. 

For both r = 2.0 and 1.4 we find best-fit mixing pa- 
rameters sin 2 29 = 1 and log(Am 2 /eV 2 ) rs -5.5. The 
absolute maximum of the likelihood is Aln(£ max ) » 3.7 
and 1.6, respectively, relative to the no- mixing case. A 
local maximum with Aln(£ max ) ss 1.4 (0.4) is found for 
sin 2 29 w 0.8 and log(Am 2 /eV 2 ) w -5. The largest in- 
crease of the maximum likelihood occurs for the largest 
relative F M temperature r = 2.0. The corresponding 
best-fit SN parameters and expected signal characteris- 
tics are listed in Tab. [| They are far away from theoret- 
ical predictions so that the apparent improvement of the 
likelihood is obtained at the price of a conflict with SN 
theory. 

Therefore, as in Sect. IV. B we next take the opposite 
point of view and assume that SN theory is roughly cor- 
rect so that we should keep Eh fixed at 3 x 10 53 erg. In 
the first analysis we allow Tv e to float for a fixed Am 2 
and sin 29 . In Fig. we show the relevant contours of 



the maximum likelihood relative to the no-mixing case. 
Again, shaded areas correspond to a diminished maxi- 
mum likelihood. As in Fig. ^ the maximum likelihood has 
an absolute maximum for t = 2.0(1.4), sin 2 29 = 1 and 
log(Am 2 /eV 2 ) w -5.7(5.6) with Aln(£ max ) « 1.4(1.1). 
A local maximum with Aln(£ max ) ps 0.8(0.3) is found 
for sin 2 29 ~ 0.8 and log(Am 2 /eV 2 ) ps -5. A similar 
effect occurred in Fig. || where the SN binding energy was 
also allowed to float. 

Next, we hold both spectral characteristics fixed, to 
wit E h = 3 x 10 53 erg and T Ve = 4.7 MeV which corre- 
sponds to the low end of the range of predicted (-E^) 
values given in Eq. ([j]). The contours of ln(£) relative 
to the no- mixing case are shown in Fig. |^ in steps of 1, 
again with t — 2.0 (upper panel) and r = 1.4 (lower 
panel). Note that all contours now represent negative 
Aln(£), i.e. diminished likelihood values. If we take the 
predicted SN parameters seriously we arrive at the same 
conclusion as in Sect. IV. B, namely that the no-mixing 
case is favored. 
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FIG. 8. Contours of hi(£) in steps of 1 relative to the 
no-mixing case. All values are negative, i.e. the maximum 
is on the line sin 2 20o = 0. The spectral parameters were 
held fixed at E h = 3 x 10 53 erg and {E- c ) = 14 MeV. The rel- 
ative Vfj, temperature r was 2.0 (upper panel) and 1.4 (lower 
panel) . 
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FIG. 10. Best-fit values for T- c and Eb and contours of 
constant likelihood which correspond to 95.4% confidence re- 
gions. In each case a joint analysis between both detectors 
was performed with sin 2 2G>o = 0.8 and Am 2 = 10~ 5 eV 2 . 
The curves are marked with the relative 77 M temperature r. 
The hatched region corresponds to the theoretical predictions 
of Eqs. (§) and (§). 



Finally, we may suppose that future experiments will 
establish the large-angle MSW solution of the solar neu- 
trino problem, i.e. that the mixing parameters lie within 
the indicated contour of Fig. |j| Specifically, we choose 
the parameters sin 2 29o = 0.8 and Am 2 = 10~ 5 eV 2 with 
t = 1.0, 1.4, 1.7, and 2.0 where r = 1.0 corresponds 
to no mixing. As in Sect. IV. B we find that the 95.4% 
confidence regions barely touch the lowest predicted v e 
energies only in the no-mixing case. However, because 
of the Earth effect the other cases yield a serious conflict 
only when the relative temperature is assumed to be 
large. 
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FIG. 9. Mixing parameters favored by the MSW solutions 
of the solar neutrino problem and those excluded by the ab- 
sence of an observed day-night effect at Kamiokande. (Con- 
tours according to Hata and Haxton (5)].) 



V. DISCUSSION AND SUMMARY 



We have studied the impact of neutrino mixing on 
the interpretation of the SN 1987A neutrino signal, fo- 
cussing on those parameter regions which are favored by 
the oscillation interpretation of the solar neutrino prob- 
lem. For these purposes the small-angle MSW solution is 
equivalent to no mixing at all because only large vacuum 
mixing angles lead to significant modifications of the an- 
tineutrino signal from a SN. In agreement with previous 
authors we find that in the no-mixing case the inferred 
neutron-star binding energy Ey> and spectral T> e tempera- 
ture are consistent with theoretical predictions, but only 
marginally so with regard to T Ue ; the 95.4% confidence 
contour in the E^-T^ plane just barely touches the pre- 
dicted range of average v e energies given in Eq. (Q) . 
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Neutrino oscillation effects lead to a partial swap of the 
F e with the stiffer spectrum. The data already point 
to lowish neutrino energies, especially at the Kamiokande 
detector, so that even a partial spectral swap aggravates 
the disagreement between the predicted and experimen- 
tally inferred neutrino energies. 

For the large-angle MSW solution the regeneration ef- 
fect in the Earth always goes in the direction of partly un- 
doing the swap caused by the adiabatic oscillation in the 
SN envelope. Therefore, the 95.4% confidence contour in 
the Sb-Tx7 e plane may be shifted only by a small amount, 
depending on the exact mixing parameters, and depend- 
ing on the relative temperature (Fig. |l0| ). Even for 
t = Tv /Tp e = 2.0 it would be difficult to claim a truly 
convincing conflict between observations and SN theory. 
Of course, the true value of r is not known. Put another 
way, if the large-angle MSW solution would be borne 
out by future solar neutrino experiments, the observed 
SN 1987A signal would have to be taken as evidence for 
a soft Vfj, spectrum relative to the v e one. 

The solar "vacuum solution" corresponds to a very 
small Am 2 for which the SN oscillations are not adia- 
batic, i.e. we have vacuum oscillations between the SN 
and here, and no regeneration effect in the Earth. In this 
case the tension between the predicted and observation- 
ally inferred SN neutrino spectra would be too significant 
to ignore, i.e. one would be forced to take the possibil- 
ity seriously that the v ^ spectra and/or V e spectra are 
softer than had been thought previously. Conversely, if 
one could show that theoretical spectral predictions were 
accurate within the claimed range of possibilities, then 
one would have to agree with the findings of Smirnov, 
Spergel, and Bahcall Q that the solar vacuum solution 
is incompatible with SN 1987A data. The conclusion of 
Kernan and Krauss [|l4| that large mixing angles were 
actually favored by the data can be upheld only if one 
ignores current theoretical predictions of the SN spec- 
tra. In this case, indeed, the likelihood function has a 
maximum for large mixing angles. 

At the present time we would argue that the theoreti- 
cal predictions of SN neutrino spectra is not well enough 
established to achieve a convincing selection between one 
of the three solutions of the solar neutrino problem. We 
note, for example, that current numerical calculations of 
the nonelectron-flavored neutrino spectra are based on 
energy-conserving neutrino-nucleon scatterings between 
their energy sphere and transport sphere in a SN core. 
However, nuclear recoils as well as inelastic modes of en- 
ergy transfer may soften these spectra in a nonneglibile 
fashion p2| . There may be other novel effects which mod- 
ify these spectra. 

Therefore, we believe that one should view the solar 
neutrino experiments as one method for shedding new 
light on SN neutrino spectra. Of course, the most inter- 
esting case would be if one of the large-angle solutions 
would obtain as they would provide nontrivial new infor- 
mation on the spectral characteristics of the SN 1987A 
neutrinos. 



NOTE ADDED 



After this paper had been submitted for publication, 
a new study has appeared where the impact of gravita- 
tional fields on the phase evolution of oscillating neutri- 
nos is investigated ^3|. We believe that for the range 
of mixing parameters and oscillation paths considered in 
our paper the gravitationally induced phases do not cause 
an observable effect. 
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